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, The remnant resulting from the merger of two neutron stars produces 

C$ ■ neutrinos in copious amounts. In this paper we present the neutrino emis- 

sion results obtained via Newtonian, high-resolution simulations of the co- 
' alescence event. These simulations use three-dimensional smoothed particle 

b ! r 

hydrodynamics together with a nuclear, temperature dependent equation of 
state and a multi-flavour neutrino leakage scheme. Wc present the details of 
our scheme, discuss the neutrino emission results from a neutron star coales- 
cence and compare them to the core-collapse supernova case where neutrino 
emission has been studied for several decades. The average neutrino energies 
are similar to those in the supernova case, but contrary to the latter, the lu- 
minosities are dominated by electron-type antineutrinos which are produced 
in the hot, neutron-rich, thick disk of the merger remnant. The cooler parts of 
this disk contain substantial fractions of heavy nuclei, which, however, do not 
influence the overall neutrino emission results significantly. Our total neutrino 
luminosities from the merger event are considerably lower than those found in 
previous investigations. This has serious consequences for the ability to pro- 
duce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted 
preferentially along the initial binary rotation axis, an event seen "pole-on" 
would appear much brighter in neutrinos than a similar event seen "edge-on" . 
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1 INTRODUCTION 



Binary pulsars such as the famous PSR 1913+16 are fascinating laboratories for extreme 
physics. Soon after its discovery it was realized that the orbit of PSR 1913+16 is decaying 
due to energy constantly leaking out of the system in the form of gravitational waves (Taylor 
1994) and therefore making the final coalescence an inescapable consequence. This merger 
event holds promises for diverse as gamma-ray bursts (Paczyhski 1986, Eichler et 

al. 1989, Narayan et al. 1992), ground-based gravitational wave detection (Abramovici et 
al. 1992, Kuroda et al. 1997, Bradaschia et al. 1990, Danzmann 1997) and the formation 
of rapid neutron capture elements (Lattimer & Schramm 1974, Lattimer & Schramm 1976, 
Symbalisty & Schramm 1982, Eichler et al. 1989, Rosswog et al. 1999, Freiburghaus et al. 
1999). 

The involved physics of the event is in almost every aspect "exotic" : the neutron star fluid 
moves in and determines the dynamical, curved space-time; in the centers of the stars and 
the resulting merger remnant the baryon densities reach multiples of the nuclear saturation 
density, p s = 2.5 • 10 14 gem -3 ; nuclear reactions proceed via extremely neutron-rich and 
short-lived isotopes; and the initial neutron star magnetic fields are expected to be ampli- 
fied during the merger to a strength, B ~ 10 17 G, so that their feed-back on the fluid flow 
becomes dynamically important (Thompson & Duncan 1993, Thompson 1994, Kluzniak & 
Ruderman 1998, Rosswog & Davies 2002). 

Neutron star mergers represent a severe challenge for computer simulations. The event is 
genuinely multidimensional and unlike, for example, with core-collapse supernovae (SNe), 
there are no basic open questions that could be first addressed in restricted dimensions, as 
e.g. the robustness of the delayed neutrino-driven supernova mechanism. Nevertheless, the 
two events have several common aspects: compact objects are formed in the center of the 
event and huge amounts of gravitational binding energy, of order 10 53 erg, are released in 
form of neutrinos making neutrino physics a key ingredient of both scenarios. The material in 
the innermost layers of both configurations is very dense, neutron rich, and neutrino opaque. 
Most neutrinos are radiated from a hot and thick accretion disk in the neutron star merger 
case, and from a shock heated mantle in the standard supernova scenario. The neutrino 
emission and absorption are the key features in the picture of a neutrino-driven supernova, 
which has been sketched in the sixties (Colgate & White 1966, Arnett 1967, Schwartz 1967), 
then refined in the mid-eighties (Wilson 1971, Arnett 1977, VanRiper & Lattimer 1981, Bow- 
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ers & Wilson 1982, Wilson 1985, Bethe & Wilson 1985, Bruenn 1985), and still continues 
to be controversially discussed and improved by many researchers. One result of this effort 
is the emergence of sophisticated neutrino transport schemes (although currently restricted 
to low spatial dimensions) to address the viability of the neutrino-driven supernova model. 
Starting with leakage schemes that considered only neutrino emission (VanRiper & Lattimer 
1981, Baron et al. 1985), multi-group flux limited diffusion approximations (Arnett 1977, 
Bowers & Wilson 1982, Bruenn 1985, Myra et al. 1987, Bruenn et al. 2001) have been devel- 
oped that take the energy spectra and a truncated expansion in the propagation direction 
between emission and absorption into account. While multidimensional simulations relying 
on transport approximations with externally imposed neutrino fluxes or spectra (Herant et 
al. 1994, Burrows et al. 1995, Janka & Miiller 1996, Mezzacappa et al. 1998, Fryer & Warren 
2002) throwed a bridge between simulation and observation, the traditional investigations 
in spherically symmetric geometry proceeded to solutions of the complete Boltzmann trans- 
port equation in stellar core collapse (Mezzacappa & Bruenn 1993) and postbounce evolution 
(Rampp & Janka 2000, Mezzacappa et al. 2001, Burrows & Thompson 2002), including full 
general relativity (Liebendorfer et al. 2001). The coalescence of neutron stars occurs on a 
much shorter time scale, of order milliseconds, compared to the shock revival in a supernova, 
which is believed to take several tenths of a second. Although weak interactions provide an 
important mechanism for the cooling of the disk that is opaque to all forms of electro- 
magnetic radiation, they do not allow for dramatic changes in the temperature and electron 
fraction, at least not on time scales accessible to current numerical simulations. Multidimen- 
sional kinematics seems to remain the dominant ingredient of neutron star mergers. 
Due to the complexity of the event simulations are still divided into two classes: either fo- 
cussing on the strong-field gravity aspect (Oohara & Nakamura 1997, Ayal et al. 2001, Faber 
& Rasio 2000, Faber et al. 2001, Wilson et al. 1996 , Baumgarte et al. 1997, Oechslin et al. 
2002, Shibata 1999, Shibata & Uryu 2000, Shibata & Uryu 2002) thereby sacrificing possibly 
important microphysics or exploring microphysics but using essentially Newtonian gravity 
(Ruffert et al. 1996, Ruffert et al. 1997, Ruffert & Janka 2001, Rosswog et al. 1999, Rosswog 
et al. 2000, Rosswog & Davies 2002, Rosswog & Ramirez- Ruiz 2002). Neutrino physics has, 
to our knowledge, so far only been included in the simulations of Ruffert et al. (see Ruffert 
& Janka 2001 and references therein) and in Rosswog & Davies (2002). 
In this paper, we detail on our neutrino leakage scheme that has been used in our high- 
resolution, three-dimensional simulations of merging neutron stars and report on the corre- 
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sponding neutrino emission results. Our leakage scheme is meant to join the current state- 
of-the-art for this specific, three-dimensional application where simplicity and numerical 
efficiency are valuable assets. It is not supposed to compare with much more elaborate 
(but low- dimensional) transport schemes necessary for quantitative statements about possi- 
bly neutrino-driven supernovae. Knowing how important the stiff energy dependence of the 
weak interactions is in the supernova, we design the leakage scheme to avoid the usage of 
mean energies for the determination of neutrino source functions or opacities. We determine 
for each neutrino energy separately a production rate and a diffusion time scale. The latter 
depends on a non-local estimate for the optical depth from which we extract the explicit 
energy dependence. The rates for the production of new neutrinos and the diffusion of neu- 
trinos from local equilibrium are then analytically integrated over energy. The smoothed 
minimum of production and diffusion rates is used as leakage source in the hydrodynamics 
equations. We apply this procedure separately for the lepton number and energy transfer. 
In Section 2 we will summarize previous results, in Section 3 we report on the neutrino 
emission results from our merger simulations. The summary and a discussion of the results 
is provided in Section 4 and the details of the neutrino treatment are given in the Appendix. 



2 BASIC MODEL FEATURES AND PREVIOUS RESULTS 

We have performed a set of high-resolution simulations of the last inspiral stages and the final 
coalescence of a double neutron star system. Large parts of the model and the hydrodynamic 
evolution have been described in detail in Rosswog & Davies (2002), hereafter referred to as 
paper I. The numerical runs analyzed in this paper are, apart from additional test runs, the 
same as those described in paper I. Here we focus on the parts of the model and the results 
that are related to the emission of neutrinos. 

Keeping in mind its decisive role for the (thermo-) dynamical evolution of the merger 
event (see e.g. Rosswog et al. 1999, Rosswog et al. 2000) we use an equation of state (EOS) 
for hot and dense nuclear matter. Our equation of state is based on the tables provided 
by Shen et al. (1998a, 1998b). We have added the lepton and photon contributions, and 
extended it smoothly to the low-density regime with a gas consisting of neutrons, alpha 
particles, electrons, positrons and photons. For details concerning the EOS we refer to paper 
I. The Newtonian self-gravity of the fluid is calculated efficiently via a binary tree (Benz et 
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al. 1990). The back- reaction forces that emerge from the emission of gravitational waves are 
added in the point-mass limit of the quadrupole approximation. 

To solve the equations of hydrodynamics for the neutron star fluid we have applied the 
smoothed particle hydrodynamics method (SPH; e.g. Benz (1990) or Monaghan (1992)). It is 
a widespread misconception that SPH is viscous "by nature" and thus necessarily introduces 
artefacts in simulations of low- viscosity flow. First, the degree of viscosity present in SPH 
is, as in every numerical scheme, a function of the numerical resolution. The components 
of the SPH artificial viscosity tensor scale to leading order proportional to the smoothing 
length h, which tends to zero with increasing resolution. The standard form of the SPH 
artificial viscosity tensor (e.g. Monaghan (1992)) is known to introduce spurious forces in 
pure shear flows. We have applied a switch suggested by Balsara (1995) which suppresses 
these forces in case of pure shear and reproduces the original form in case of shocks. A 
further improvement concerns the artificial viscosity parameters, usually called a and (3: they 
are made time dependent (as suggested in Morris and Monaghan 1997) and an additional 
differential equation is solved to determine their values. In the absence of shocks these 
values are negligible, if a shock is detected the parameters rise to their standard values. This 
artificial viscosity treatment is described and tested in detail in Rosswog et al. (2000). 
To quantify the amount of viscosity in our current simulations we have estimated the effective 
a-viscosity present in the disk of the merger remnant. The effective a- viscosity is ass cc h/H, 
where H is the thickness of the disk, and therefore depends on how well-resolved the vertical 
disk structure is. We found very low numerical values, ass ~ 10~ 3 for the disks in our models 
and even lower values in the better resolved central regions of the remnant. 
The whole code is parallelized for shared-memory architecture and obtains an excellent 
speed-up for up to ~ 100 processors. In a typical application with several 10 5 particles a 
speed-up of 55 is obtained on 60 processors. 

We follow the system evolution from an initial separation of ~ 3-R ns , where R QS is the 
radius of an isolated neutron star, for approximately 15 ms. From the chosen initial sepa- 
ration it takes the neutron stars only a few milliseconds to merge. They leave behind an 
extremely massive central neutron star (~ 2.4 M ), surrounded by a hot and dense, shock- 
heated inner disk region (with temperatures T ~ 3 MeV, densities p ~ 10 12 gem -3 and a 
mass Mdisk ~ 0.2 M ) and rapidly expanding debris material. 

The central neutron star is strongly differentially rotating, most pronounced in the generic 
case without initial spin. Since differential rotation allows the central parts to spin extremely 
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fast without the (slower rotating) outer parts of the object reaching the mass shedding limit, 
a substantially higher maximum mass can be stabilized. A recent investigation (Lydford et 
al. 2002) using polytropic equations of state finds values of (M r — M nr )/M nr up to 1.8 for 
soft EOSs and for the polytrope closest to our nuclear EOS they find (M r — M nr )/M nr up 
to 0.6 which corresponds to masses well beyond the total binary mass of 2.8 M (M r is 
the maximum mass for a differentially and M nr . the maximum mass of a non-rotating star). 
We therefore expect an extremely massive, hot neutron-star-like object to be formed in the 
center of the merger remnant whose lifetime is determined by the time it takes to get rid of 
the rotational support. Although a conclusive answer to this point cannot be given from the 
current calculations (since they are essentially Newtonian and some of the physics ingredi- 
ents like the high-density part of the equation of state are to date only poorly known), we 
estimate that the neutron star might remain stable for many dynamical time scales. This 
time scale may be long enough to allow the magnetic seed fields to be amplified to enormous 
field strengths (~ 10 17 G; Thompson and Duncan 1993). If one assumes magnetic dipole 
radiation to drive the system towards black hole formation, even time scales of months can 
be easily obtained without stretching the involved parameters beyond reasonable limits. The 
exact time scale between the merger event and the (probable) final black hole formation may 
depend quite sensitively on the details of the specific merger event. For a further discussion 
see Rosswog & Davies (2002). 

The debris around the central object exhibits an interesting flow-pattern: material that has 
previously been centrifugally launched into eccentric orbits, and thereby cooled by expansion 
and neutrino-emission, is returning towards the central object. This cool (T < 0.5 MeV), 
equatorial inflow produces a butterfly-shaped shock front when it encounters material that is 
still being shed from the central object. In this way a hot flow is driven in vertical direction. 
The resulting disk is very thick with a height comparable to its radial extension. 

In the present paper we will report on these simulations with a focus on the neutrino 
emission that goes along with a neutron star coalescence. 

3 NEUTRINO EMISSION FROM NEUTRON STAR MERGERS 

Under the conditions of a neutron star merger neutrinos are produced copiously and they 
provide the most efficient cooling mechanism for the dense, shock- and shear- heated neu- 
tron star debris. In addition, the related weak interactions determine the compositional 
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Table 1. Summary of the different runs, ao : initial separation; v : neutrino physics; T s ; m : simulated duration; M1/M2 : 
masses in solar units; # part.: total particle number 



run 


spin 


Mi 


M 2 


# part. 


a () [km] 


V 


T sim [ms] 


remark 


A 


corot. 


1.4 


1.4 


207,918 


48 


no 


10.7 




B 


corot. 


1.4 


1.4 


1,005,582 


48 


no 


10.8 




C 


irrot. 


1.4 


1.4 


383,470 


48 


yes 


18.3 


gen. case 


D 


corot. 


1.4 


1.4 


207,918 


48 


yes 


20.2 


lower i/-limit 


E 


irrot. 


2.0 


2.0 


750,000 


48 


yes 


12.2 


upper ^-limit 


F 


corot. 


1.4 


1.4 


20,886 


52.5 


yes 


15.1 


spur, ^-emission ? 



evolution via the electron fraction Y e that is altered by charged- current reactions such as 
electron and positron captures. The enormously temperature dependent weak interaction 
processes can exhibit in some parts of the flow very short time scales, \Y e /Y e \ ~ 1CT 6 s, 
which is well below the dynamical time scale of a neutron star, r dyn = (Gp) -1 / 2 ps 2 • 1CT 4 
s, while they are essentially infinite in other parts of the flow. Therefore neither the as- 
sumption of an instantanous beta equilibrium nor frozen Y e values are justified. Since in the 
dense parts of the hot, merged configuration the neutrino mean free paths are of the order 
A ~ 0.75 m (5 • 10 14 gem" - 3 /p) (10 MeV/T) , where p is the matter density and T is the 
temperature, the interaction of the neutrinos with the ambient matter has to be accounted 
for. Here and in the rest of the paper we measure temperatures in energy units, i.e. ks = 1. 
A full Boltzmann neutrino transport in the context of the three-dimensional modelling of 
the event is beyond the current state-of-the-art and computational resources. But since the 
simulated physical time scales are of the order of 10 milliseconds and neutrino momentum 
transfer is expected to be unimportant we consider a detailed neutrino leakage to be an 
important step towards reliable physical models of the event. 

We consider three neutrino flavours: electron neutrinos, v e , electron anti-neutrinos, P e , and 
the heavy-lepton neutrinos, v^v^v T ^v T ^ which are collectively referred to as v x . The basic 
idea of our leakage scheme is to provide a physical limit via diffusion rates. This guarantees 
the limitation of the neutrino production to the amount that is able to stream away. In the 
opaque regime the neutrinos therefore escape only on a diffusion time scale and in the trans- 
parent regions they leave their production site essentially without any further interaction 
with the surrounding matter. 

The dominant neutrino processes in our context are the charged-current lepton capture 
reactions on nucleons, electron capture (EC) 

e~+p-^n + v e (1) 
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and positron capture (PC) 

e + + n -> p + P e (2) 

which produce electron flavour neutrinos and the "thermal" , pair producing reactions, pair 
annihilation 

e~ + e + — > i/j + Pj (3) 
and plasmon decay 

7 -> z/i + Pi, (4) 

which produce neutrinos and anti- neutrinos of all flavours, i/j and Pj. The latter process 
dominates in the strongly electron-degenerate regime. We disregard electron captures onto 
nuclei since these reactions would require detailed information about the nuclear shell struc- 
ture which is not available for these nuclei. But these captures are not expected to be 
important in our case, since the regimes where the dominant neutrino emission takes place 
are almost completely photodisintegrated (see below). We further neglect neutral-current 
nucleon-nucleon bremsstrahlung as neutrino production process. This process has recently 
received attention in the supernova context (Thompson, Burrows and Horvath 2000). It may 
be possible that this process is locally important, but recent investigations (Keil et al. 2002) 
including this and other reactions for the supernova case only found overall changes of the 
order 10 % . Since accurate emission rates are difficult to obtain (due to the poor knowledge 
of the nucleon-nucleon potential and due to uncertainties in the magnitude of many-body 
effects) and we do not expect effects larger than the uncertainties inherent in our leakage 
scheme, we decided to ignore this process. 

To determine the number and energy diffusion rates based on neutrino opacities we take 
into account the scattering off nucleons, 

Vi + {n,p} -> Vi + {n,p}, (5) 
coherent neutrino nucleus scattering, 

Vi + A -> Vi + A, (6) 
and neutrino absorption by free nucleons, 

v e + n — > p + e~ (7) 
v e +p^n + e + . (8) 



High Resolution Calculations of Merging Neutron Stars II: Neutrino Emission 9 

For the details of the implementation of the reactions and the leakage scheme we refer to 
the Appendix. 

3.1 Total luminosities, mean energies 

Neutron stars are expected to heat up during inspiral by tidal interaction to temperatures 
of the order 10 8 K (Lai 1994). At these temperatures no significant neutrino emission will 
occur. In order to test for the amount of "spurious" emission of neutrinos in our simulation 
due to the unavoidable numerical heat up of the completely degenerate stars, we perform 
a test run (a listing of the different runs is provided in Table HJ). We prepare a corotating 
equilibrium binary configuration just outside the last stable orbit by relaxing the two neutron 
stars in their mutual gravitational field. Subsequently we follow their dynamical evolution for 
approximately 50 neutron star dynamical time scales while they revolve on perfectly circular 
orbits around their common center of mass. Only self-gravity and hydrodynamic forces are 
considered, no initial radial velocities are applied and the gravitational wave backreaction 
forces are switched off. We find that the total neutrino luminosity reaches a stationary level 
of ~ 3- 10 49 erg s -1 , see Fig.[Q which is four orders of magnitude below the peak luminosities 
of the full merger calculation and therefore completely negligible. 

The overall neutrino emission properties of the full merger calculations are shown in Figs. 
El to The total neutrino luminosities (left panels) are calculated summing up all particle 
contributions and the rms energies for each neutrino flavour are calculated according to eq. 
(jA4|) . For the corotating case, run D, substantial neutrino emission sets in later than in the 
cases without initial spin. The explanation for this is twofold: on the one hand this run 
starts out from numerically exact initial conditions while the non-rotating cases, run C and 
E, suffer an accelerated inspiral due to the start with initially spherical stars. On the other 
hand, neutrino emission only becomes important once the thick torus around the central 
high-density part has formed. Again this process takes longer for the corotating case since 
(due to the larger initial angular momentum) the torus-forming matter is initially launched 
into wider orbits. The average neutrino energies reach peak values soon after the stars have 
first come into contact. The reason for this is that the neutrinos from the hot debris material 
can, at this stage, escape without having to pass through any optically thick matter. 
The total neutrino emission seems to have reached a roughly stationary level (except for 
maybe run C) by the end of the simulation. The total luminosities range from ~ 10 53 erg/s 
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for the smoothly merging corotating case over ~ 2 ■ 10 53 erg/s for the irrotational case with 
twice 1.4 M to ~ 4 • 10 53 erg/s for our extreme case with 2 x 2.0 M Q and no initial spins. 
We regard run D as a lower limit which is unlikely to occur in nature (Bildsten & Cutler 
1992, Kochanek 1992) and run C as the generic case since the observed neutron star binary 
systems have masses close to 1.4 M (Thorsett & Chakrabarti 1999) and are expected to 
have a very slow individual spin at the merger stage. The extreme case run E has been 
performed in order to explore the upper limit on the neutrino emission from the merger 
event. 

The mean energies are ~ 8 MeV for the electron-type neutrinos. This is below the typical 
values found in core collapse supernovae. The material in the nascent protoneutron star has 
first to deleptonize on a neutrino diffusion time scale before the electron fractions are as 
low as in the neutron star merger event, where beta-equilibrium has been established in 
the individual neutron stars long before the coalescence. Hence, the material in the super- 
nova is more electron degenerate at comparable densities. Electrons are therefore captured 
from higher Fermi energies and produce electron neutrinos with a harder spectrum in the 
supernova case. The situation is different for the electron antineutrinos. We find rms ener- 
gies around ~ 15 MeV, quite comparable to rms energies in the supernova case. The lower 
electron degeneracy in the neutron star merger favours the population of positrons, whose 
chemical potential has to balance the electron chemical potential because of pair equilib- 
rium. The high positron abundance in combination with the neutron rich matter leads to 
more positron capture events on free neutrons than in the supernova. This results in higher 
electron anti-neutrino luminosities. The electron anti-neutrino luminosity from the remnant 
reaches up to ~ 1.5 • 10 53 erg/s, it provides the main cooling mechanism of the hot accretion 
disk. The heavy lepton neutrinos reach rms energies of ~ 20 to ~ 25 MeV. This is compa- 
rable to the supernova rms energies. Their luminosity, however, tends to be smaller in the 
neutron star merger case because the temperatures in the high density regimes, where the 
heavy lepton neutrinos emerge, are manifestly lower (see below). 

To characterize the physical conditions of the emission region of each neutrino flavour we 
calculate average quantities, X u ., weighted by the z/j-number production rates per particle, 
Rl{j, where X stands for p, T, Y e and /i e , given by 



k = m- f 3 . (9) 



These quantities are displayed in Table El Electron neutrinos and anti- neutrinos are emitted 
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under similar conditions, typically at densities around 10 12,5 gem -3 , temperatures of 4-5 
MeV (anti-neutrinos at slightly higher values) and a Y e below 0.1. The heavy lepton neutrinos 
are emitted at substantially higher densities (log(p) ~ 13) and temperatures (T « 9 MeV). 
Figure 03 shows a comparison between the conditions encountered in neutron star mergers 
and those of SNe. In the upper panel we show SPH-particle densities and temperatures (black 
dots, every 20th particle is shown) for our generic run C. Particles with peak luminosities 
are indicated with special symbols (these peak values are not to be confused with the average 
properties mentioned above). Filled circles indicate particles which emit u e at a luminosity 
in excess of 10 % of the maximum particle z/ e -luminosity, squares mark the corresponding 
particles for v e emission and triangles refer to v x . Due to the steeper temperature dependence 
{Qv x oc T 9 ) a lower threshold (3.5 %) has been chosen for the v x in order to display roughly 
the same number of particles. The corresponding plot for Y e as a function of log(p) is shown 
in the second panel of Figure El We compare the state of these fluid elements to the state of 
fluid elements in a simulated postbounce evolution of the core of a 13 M progenitor star. 
Due to the limitation to spherical symmetry in this simulation with Boltzmann neutrino 
transport (Liebendorfer et al., 2002), the fluid elements form a solid line. In the supernova 
case, at 100 ms after bounce, we find the peak emission at densities of 10 106 , 10 11 ' 3 , and 
10 12 ' 6 g/cm 3 respectively. This is in agreement with the high emission regions identified in 
the neutron star merger for the heavy lepton neutrinos and the electron anti-neutrinos. The 
peak emission of electron neutrinos in the neutron star merger appears to occur at slightly 
higher densities (~ 10 11 g/cm 3 ). We attribute this to the less pronounced compression and 
deleptonization of infalling matter at low densities in the rotating accretion disk if compared 
to the failed explosion of a non-rotational supernova simulation. The heavy lepton neutrinos 
stem in both SN and neutron star merger from similar densities. 

To analyze the importance of the e + /e~-capture reactions, eqs. (0) and (J2J), versus the 
pair producing reactions eqs. (jHJ) and we perform a post-processing experiment. We 
take one time-slice of our generic run, run C, at t= 14.1 ms and use the pair and plasma 
neutrino reactions as the only emission processes (i.e. the capture reactions are artificially 
switched off). In this case the luminosity in electron-type (anti-) neutrinos is only ~ 10% 
of the previous values, indicating that a major contribution stems from the lepton capture 
reactions. 
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Table 2. Typical properties of emission region, densities are in gem -3 , temperatures and chemical potentials in MeV. 



run 


log{p) Vc 
























C 


12.6 


12.6 


13.2 


4.2 


5.5 


8.9 


0.072 


0.072 


0.13 


12.1 


13.2 


23.7 


D 


12.4 


12.2 


13.5 


4.1 


5.0 


6.6 


0.083 


0.095 


0.10 


11.6 


13.0 


35.1 


E 


12.1 


12.4 


12.9 


5.0 


5.7 


9.1 


0.140 


0.085 


0.12 


9.6 


13.0 


22.1 



3.2 Emission geometry: disk versus central object 

In Fig. |U1 we plot the neutrino energy (sum of all flavours) per time and volume for run C 
(upper two panels), run D (intermediate two panels) and run E (the two lower panels). The 
left column of panels shows the emission in the orbital plane, while the vertical emission ge- 
ometry (azimuthally averaged) is displayed in the right column. Note that the emission per 
time and volume from the hot, but extremely dense central objects is completely negligible, 
roughly two orders of magnitude lower than that coming from the most luminous parts of 
the disk. In paper I we had mentioned the butterfly-shaped temperature distribution in the 
XZ-plane that results from cool inflow being shock-heated hitting the inner parts of the disk, 
see Figure 15 in paper I. This pattern is also reflected in the neutrino emission geometry, 
see right column in Fig. |H1 

3.3 Opacity sources: Importance of heavy nuclei 

Our scheme accounts for the coherent scattering of neutrinos off heavy nuclei, for details 
we refer to the Appendix. We had realized that, despite the high temperatures encountered 
in the disk, matter finds it energetically favorable to form a non-negligible mass fraction 
of heavy nuclei (paper I). Due to the (approximate) proportionality of the scattering cross- 
section to the square of the nucleon number of the heavy nucleus, see (jA18|) . nuclei could 
possibly dominate as an opacity source. To estimate how important the heavy nuclei really 
are for the neutrino emission, we perform the following test. We take one time slice (t= 
14.1 ms) of our generic case, run C, update the neutrino grid (see Appendix) and then 
calculate with these opacities the properties of the emitted neutrinos. In one case we use 
-like in the dynamical simulation- the full set of abundances given by the EOS for both the 
emission and absorption/scattering processes and in the other case we assume the matter 
to be completely dissociated into nucleons, i.e. the mass fractions are given by 

x p = Y e , x n = 1 - Y e , x a = 0, x h = 0. (10) 
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We find almost exactly the same numbers for both the mean energies and the total lumi- 
nosities, maximum deviations in the (more sensitive) total luminosities are below 5 %. 
The reason for this lies in the geometry of the heavy nucleus distribution. In Fig. Owe show 
the azimuthally averaged values of the heavy nucleus mass fraction, Xh, of runs C, D and 
E; these values are shown for matter with densities above 10 10 gem -3 , below that density 
matter is transparent to neutrinos. Nuclei are present in the cool, equatorial inflow regions 
identified in paper I, see Figure 15 in Rosswog and Davies (2002). The butterfly-shaped 
temperature distribution is also reflected in the nucleus mass fraction. It is interesting to 
note that despite the extreme temperatures in the central object a thin, nuclear crust can 
survive in our coolest case (run D). The hottest case (run E) is essentially free of heavy 
nuclei. By comparing Fig.[7|with Fig. right column, it becomes obvious that the neutrinos 
from the most luminous regimes can escape in each case vertically without having to pass 
through material containing an interesting amount of heavy nuclei. Therefore the influence 
of the heavy nuclei onto the total luminosity and the mean energies is negligible. 

3.4 Optical depths, neutrino-" spheres" 

To illustrate how the opaque matter is distributed in the merger remnant we plot in Figures 
E] to HOI contours of the spectrally averaged neutrino optical depth (see eq. (1A23J) and QASjl ). 



for all neutrino species at the end of run C, D and E. The F n are the standard Fermi integrals, 
see Appendix A. Consistent with our neutrino treatment we have used the equilibrium values 
(see eq. (|A5|i ) for the degeneracy, T], of the v e jv e and have assumed a vanishing degeneracy 
parameter for the v x . In each of the Figures the uppermost panel shows contours of the optical 
depth of the electron neutrinos, the middle panels refer to the electron anti-neutrinos and 
the remaining ones to the heavy lepton neutrinos. 

The debris matter is most opaque to the electron-type neutrinos which in addition to 
scattering are also absorbed onto the copiously available free neutrons. Electron-type anti- 
neutrinos see matter less opaque and matter is most transparent to the v x . We also show (as 
the thick line) the "neutrino-sphere", defined as the locus with r Vi = 2/3. For the v e and the 
v x the neutrino-spheres almost coincide, with radial extensions of ~ 70 km and peak heights 
of ~ 20 km, since both neutrino types suffer essentially the same interactions (scattering 
events; in the extremely neutron-rich debris absorption onto free protons is only a minor 




(11) 
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correction). Due to their additional opacity sources the v e decouple substantially further 
out, at radial extensions of ~ 105 km with peak heights of ~ 35 km. 

The optical depths at height z= 0, i.e. in the orbital plane, are shown in Figure ITT1 (from top 
to bottom: v ei v e and v x ). In the central object values up to several times 10 4 are reached, 
beyond ~ 130 km matter is essentially transparent to neutrinos of all types, i.e. r u . < 0.1. 
It is interesting to note that it is only in the central object {r Vi > 10) that neutrinos are 
really trapped, see Figures IH1 to ITTI At the edge of the central object, at distances of ~ 30 
km from the origin, the optical depth drops rapidly, but then only decreases very slowly 
throughout the disk (~ 30 km to ~ 100 km). The whole hot torus-region is therefore in the 
semi-transparent regime. 



3.5 Directional dependence of neutrino emission 

It is consistent with our approach from eqs. (jAl|) and (jA2|) to think of the neutrinos emit- 
ted from an SPH-particle to be composed of "free neutrinos" and "diffusive neutrinos". In 
analogy to general diffusion equations we assume that the diffusive neutrino component is 
emitted in the direction of the local, negative density gradient, h = — (Vp)/| — (Vp)|. We 
use the SPH-prescription to determine this density gradient at the position of particle i: 

V Pi = Y t m j ViW ij , (12) 

3 

where rrij is the particle mass, Wij = W( ^V~ Xj ) the standard SPH-kernel (e.g. Monaghan 
1992) and hij is the arithmetic mean of the involved smoothing lengths. The free component, 
in contrast, will emit isotropically. The fraction with which the both components contribute 
to the neutrino luminosity of particle j is given by 



fdif _ ^UjJ J rloc _ V^J /io\ 

JvuQi ~ n dif dI1Cl Jvi,Qj ~ n ioc ■ \ L<i ) 



o ef o ef 

%f and fl oc Q = ^ 

n dif Jfi,Qj r\l 

V^J ^Vi,3 

It can be easily checked that fJ?Q. + /^°q — 1 by using eqs. ()A1|) and (jA2j) and that 
the fractions approach their obvious limits in the high and low-density regimes. Similarly, 
fractions of the emitted neutrino number, f v *f R . and fjf^R. can be defined as above, but with 
number emission rates per volume, R Vi ,j, rather than with energy emission rates per volume, 
Q Vu j. With these definitions the neutrino luminosity per solid angle, composed of a diffusive 
and a free component, is determined by 

A "*W- An ~27rsin(^)A^ + 4vr ' 1 } 
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Here Q Ui ,k is the z/j energy emission rate of particle k, from either "diffusive" neutrinos, 
Qtlk = fu/,Q k ■ Qt{,k or " free " neutrinos, Q l °° k = f l ° % c Qk ■ Q e J ik . In the above equation the 
j-sum extends over all the particles (since their "free" neutrinos radiate isotropically), the 
fc-sum, however, only extends over those particles that radiate into a ring of width At? in 
the ^-direction, d — A$/2 < < i? + A$/2, where is given by cos($k) = ■ e z . An 
observer that sees the merger from an angle $ with respect to the initial binary rotation axis 
(= z-axis) would thus infer an apparent luminosity of 

ZJ%>(4) = 4x1^(0). (15) 

The quantity A^fi) for our generic case, run C, is shown in Fig. ^] (the other runs yield 
similar results). The luminosity per solid angle is peaked towards the z-axis: a system ob- 
served "pole-on" (i? f& 0°) will yield a total neutrino energy flux, given by = Al ^2 , 
R being the distance to the source, that is around 20 times larger than that of a system 
that is observed "edge-on" (i9 w 90°). This preferential emission is visible for all neutrino 
flavours, but most pronounced in the case of the heavy lepton neutrinos, v x . The latter ones 
are produced in the enormously temperature dependent reactions © and (HJ) and therefore 
emerge from the hottest parts of the remnant, i.e. they are generated within or close to the 
flattened central object, where the density gradients point along the z-direction. 
To infer the -^-dependence of the average neutrino energies we use the quantitity 



(16) 



Here, the j-sum extends over all the particles; to count only the contributing particles in 
the diffusive part, the function 

2/(mn(0)A0) for M/2 <•&;<■& + &&/2 

else 

has been introduced. The distribution of the average energies, see second panel Fig. ^1 is 
relatively flat for all neutrino types. 



9,(0) 



4 SUMMARY AND DISCUSSION 

We have presented the neutrino emission results from our high-resolution simulations of the 
coalescence of two neutron stars. We find typically total neutrino luminosities of ~ 2 ■ 10 53 
erg/s with rms energies of ~ 8 MeV for electron type neutrinos, ~ 15 MeV for their anti- 
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particles and ~ 20 — 25 MeV for the \x- and r-neutrinos and their antiparticles. We have 
performed two runs that are intended to give an upper and a lower limit to the total neutrino 
luminosities: in the case of initial corotation the stars merge extremely smoothly. This goes 
along with moderately high matter temperatures and neutrino luminosities lower by a factor 
of two than in our standard case. In the other extreme case we consider the coalescence of 
two 2.0 M Q neutron stars without initial spins. This case leads to the hottest merger rem- 
nant and, correspondingly, to the highest neutrino luminosities, around 4 • 10 53 erg/s. 
The contributions of the extremely hot, but also neutrino-opaque central objects are marginal, 
typically only a few percent. Most neutrinos are produced in the debris torus around the 
central object, which exhibits temperatures well above the positron production threshold 
and which is very neutron rich (Y e ~ 0.1). These conditions favour positron captures on free 
neutrons over electron captures and therefore yield neutrino luminosities which are clearly 
dominated by the v e . The heavy lepton neutrinos contribute only ~ 10% to the total lu- 
minosity. The v e and v e are predominantly produced in electron and positron captures on 
free nucleons, only ~ 10% come from the pair and plasma process. The whole disk, with 
distances of ~ 30 km to ~ 100 km is semi-tranparent to the neutrinos; it is only within the 
high-density central object that they are really trapped (t u . > 10). 

We find qualitative agreement with the results described in Ruffert and Janka (2001) as far 
as the dominance of the v e and the hierarchies in the rms neutrino energies (e Ve < ep e < e Vx ) 
are concerned. Our total luminosities and rms energies, however, are lower than those found 
in their models, typically by a factor of ~ 2 in the luminosities and around ~ 20% in the 
mean energies. This comes in part from the fact that not exactly the same initial conditions 
are used: the 'standard' mass they use is 1.6 M Q rather than our value, 1.4 M . This mass 
difference is expected to lead to slightly increased luminosities. Another difference is the 
EOS. The EOS of Shen et al. (1998) that we use is in the density regime of 12 < log(p) 
< 14, where large fractions of the neutrino emission stems from, substantially stiffer than 
the Lattimer-Swesty EOS that Ruffert & Janka use (compare Fig. 2 in Rosswog & Davies 
2002). This leads to a less compact configuration with lower temperatures and correspond- 
ingly lower neutrino luminosities in our case. Further possibilities include a different amount 
of numerical viscosity (see Rosswog & Davies 2002 for a discussion) in both codes and maybe 
the interaction with the background medium in the simulations of Ruffert & Janka (2001; 
see their paper for a discussion of this point). Finally, the lower luminosities may also come 
from differences in the leakage prescriptions. However, as shown in the appendix, if at all, 
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our scheme tends to overestimate the luminosities. Therefore the true neutrino luminosities 
could be even lower, a fact that has serious consequences for the ability of neutron star 
mergers to produce a gamma-ray burst fireball via neutrino annihilation. This is discussed 
further in Rosswog & Ramirez- Ruiz (2002) and Rosswog et al. (2003). 
Despite the high temperatures we find areas in the disk that contain a substantial mass frac- 
tion of heavy nuclei. One might expect this to influence the neutrino luminosities, since the 
coherent scattering cross sections are oc A 2 , where A is the nucleon number of the nucleus. 
This, however, is not the case. In each of the investigated cases the most neutrino-luminous 
parts of the remnant are essentially free of heavy nuclei and the neutrinos can always escape 
via almost completely photo-dissociated matter. 

We find that the neutrino emission per solid angle is focussed towards the initial binary 
rotation axis. A merger remnant observed "pole-on" has an apparent neutrino luminosity 
that is about 20 times larger than a remnant seen "edge-on". 

A typical neutron star merger produces mean neutrino energies very similar to those result- 
ing from the core collapse of a massive star. A distinctive signature between both events is 
the strong dominance of the electron anti-neutrinos over electron neutrinos in the merger 
case. The most unique proof, however, for neutrinos coming from a neutron star coalescence 
rather than from a SN would be the nearly coincident detection of a binary "chirp" -signal 
in gravitational waves. The peak luminosity in neutrinos will be reached about 15 ms after 
the peak in the gravitational wave luminosity. 
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APPENDIX A: NEUTRINO TREATMENT 

The rates that we use in the simulations are smooth interpolations between diffusion and 
local production rates. If we denote for a given neutrino species z/j the number emission rates 
by R Vi per volume and energy emission rates per volume by Q Ui , our prescription for the 
effective rates reads 



(tdIoc \ _i 
1+ w) (A1) 



« = « C (l + ||) • (A2) 

This ansatz is similar to the one used in Ruffert et al. (1996). Here the quantities with 
the superscript "/oc" denote the locally produced rates of number and energy while the 
superscript u dif" refers to the diffusion rates that are further specified below. In the trans- 
parent regime, where the diffusion time scale T^ l J is short, and therefore Rf/ 3> R l °° and 
Qt/ ^ Q l °° an the locally produced neutrinos stream out freely. In the very opaque regime, 
where T„ % * is large, the neutrinos leak out on the diffusion time scale. Therefore both limits 
are treated correctly, the regime inbetween these limits is handled via interpolation. 
The mean neutrino energy of each SPH-particle (particle index suppressed) is then found 
from 

f = Lr^u r (A } 

* D e f V ' 

where r labels all reactions producing neutrinos of type z/j. Note, that these (mean) energies 

are used exclusively for book-keeping purposes, in all places where a dependence on neutrino 

energies occurs, we integrate cross-sections over a Fermi-distribution (see below). 

To characterize the average neutrino energies of the total system we use rms energies given 

by 



e,„ = 



^^'^fpdL, (A4) 



where j labels the SPH-particles and R e v . j is the rate of neutrino number emission of particle 
j (not to be confused with the rate per volume, R e J u j)- 
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We have tested this scheme in spherical symmetry against stationary state Boltzmann trans- 
port (Mezzacappa & Messer 1999). To this end we determined the neutrino properties for 
a frozen matter background. The background properties (p, T and Y e ) were either taken 
from neutron star merger (Rosswog & Davies 2002) or core collapse supernova simulations 
(Liebendorfer et al. 2002). While the rms neutrino energies agree within 20 % the accuracy 
of the luminosities depends on the importance of the semi-transparent regime where the 
interpolation (eqs. (jAl|) and (|A2|1 ) is applied. In the worst case we found that our scheme 
overestimates the luminosities by a factor 3-4. 



Al Free Emission Rates 

In the following we will neglect the electron mass and the nucleon mass difference, Q = m n — 
m p = 1.2935 MeV, in all the cross sections (Tubbs & Schramm 1975). This is appropriate 
for our purposes and largely simplifies the involved rate expressions. We further assume the 
neutrino temperature to be identical to the local matter temperature and, where necessary, 
we assume the neutrinos to follow a Fermi-distribution. The chemical potentials of the v x 
are generally assumed to vanish, for v e and v e we apply the equilibrium values 

/V = -A^ e = fie - fi - Q, ( A5 ) 

wherever they occur in the sequel. Here jx e is the electron chemical potential (with rest mass) 
and fi is the difference in the neutron and proton chemical potentials (without rest mass). 
Degeneracy parameters fii/T are denoted by % temperatures are always in units of energies. 
With these approximations and ignoring momentum transfer to the nucleon (Bruenn 1985) 
the electron capture rate per volume reads 

R E c = f3r ] pnT 5 F A {r ]e ),, (A6) 
with 

(3— 71 1 + ^ a2 (T (A7) 
h 3 c 2 (m e c 2 ) 2 

Here h is Planck's constant and c the speed of light, a ~ 1.25, m e is the electron mass, 
o"o ~ 1.76 • 10~ 44 cm 2 . F n is a Fermi integral given by 

Fn {z) = / (A8) 
v ; Jo e x - z + 1 ( J 

and can be efficiently evaluated via series expansions (Takahashi et al. 1978). The factor r) pn 

given by 
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^ pn exp(jl/T) — V ^ 
takes into account the nucleon final state blocking and reduces in the non-degenerate limit 
to the proton number density n p , n n refers to the neutron number density. Following the 
analogous procedure one finds for the energy emission rate 

QEc = (3v P nT 6 F 5 ( Ve ), (A10) 
and for the mean energy of the emitted neutrinos 

IW \ _ Qec _ T F 5(Ve) . 

{e ^ e °-r^- t fS)- (A11) 

The corresponding rates for positron captures read 

Rpc = PVn P T 5 F i {-r ]e ), (A12) 

Qpc = PVn P T 6 F 5 (- Ve ), (A13) 

(E^) PC = T§^\, (A14) 

where r] np is obtained from r] pn by interchanging the neutron and proton properties. 
The "thermal" processes are taken into account via fit formulae. For the energy emission 
from the pair process we use the prescription of Itoh et al. (1996). The number emission rate 
is obtained by deviding by the mean energy per neutrino pair (Cooperstein et al. 1986) 

(E^Ur = T(^ + i^\ ). (A15) 

For the plasmon decay we use the formulae of Haft et al. (1994) with 

<Mi = tU + j^A , (A16) 



where 7 = 7o-i/7r 2 /3 + r\\ and 70 = 5.565 ■ 10" 



A2 Diffusive Emission Rates 



In order to evaluate the opacities along given directions we map the particle properties den- 
sity, temperature and electron fraction on an aequidistant, cylindrical grid with coordinates 



(R, Z), where R = \/x 2 + y 2 , see Figure ITBI The assumption of rotational symmetry around 
the binary rotation axis is an excellent approximation since the main neutrino emitting re- 
gion is the hot, neutron star matter debris torus that forms around the merged central object 
(see Fig. 14 in paper I). By evaluating the EOS at each grid point the matter properties 
(like the local composition) are known and we can therefore assign a variable ( Ui (see eq. 
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()A21[) ). containing compositional information, to each grid point. The neutrino grid does not 
have to be updated at every hydro time step. We chose to update it after a small fraction 
(1/8) of the neutron star dynamical time scale, Td yn — (Gp)^ 1 ^ 2 ~ 2 ■ 10~ 4 s, which is a tiny 
fraction of the timescale on which typical disk properties change. Once all the properties on 
the grid are known, the desired values at the SPH-particle positions are found by trilinear 
interpolation. We use 400 points in radial direction and 300 points in positive Z-direction 
(symmetry with respect to the orbital plane is an excellent approximation for the systems 
under investigation). 

The dominant sources of opacity are 

(i) neutrino nucleon scattering: 

Vi + {n, p} -> v>i + {n, p} (A17) 

with a uunuc = |cr Q (^2) (Shapiro & Teukolsky 1983) and 

(ii) coherent neutrino nucleus scattering: 

v>i + A -> vi + A (A18) 

with o~ Uij A = jq°~o (J^s) A 2 (l — Z/A) 2 (Shapiro & Teukolsky 1983; sin 2 8w has been approx- 
imated by 0.25). Here A and Z are the nucleon and proton number of the average nucleus 
whose properties are stored in our EOS-table. Due to the /^-dependence of the cross sec- 
tion this process will dominate as soon as a substantial fraction of heavy nuclei is present 
(remember that the nucleon numbers in these nuclei reach values of up to ~ 400 Shen et al. 
1998). 

Electron type neutrinos additionally undergo 

(iii) neutrino absorption: 

v e + n -> p + e~ (A19) 
v e + p -> n + e + (A20) 

with a Ve>n = iif^ao {-^f (1 - f e -), where (1 - f e -) w (exp^ - F 5 (^J/F 4 (^J) + l)" 1 
and a PeiP = ±±f^a {-£^f (1 - / e+ >, (1 - f e+ ) « (exp(-r ?e - F 5 (^J/F 4 (^J) + l))" 1 . 

The local mean free path is given by (where for simplicity the spatial dependence is 
suppressed) 

K(E) = fenra r (£)j = (E%r\ (A21) 
where the n r denote the target number densities, the index r runs over the reactions given 
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above with cross-sections a T and E is the neutrino energy. The dependence of the cross- 
sections on the squared neutrino energies has been separated out in the definition of £ Ui . The 
optical depth, r, along a specified direction is then given as 

^-Ltj® (A22) 

The optical depths are evaluated along three directions from each grid point: in Z-direction 
(t*.), i.e. parallel to the rotational axis, along the outgoing diagonal (r£) and along the 
ingoing diagonal (r^.), see Fig. [T31 The finally used optical depth, (r Ui ), is the minimum of 
the three, t Vi = min(r* , t%., t% ). The quantities that are actually stored for each grid point 
j are 

XL = I U*)dx, (A23) 

where d denotes the direction and J d • dx is the integration from grid point j along direction d. 
Note that the quantity x is independent of the neutrino energy and the (energy dependent) 
optical depth is given by 

r Ui {E) = E 2 mm d (xU) = E\^, (A24) 

The diffusion rate depends on the optical depth t Vj . We base our estimates on a very 
simple, one-dimensional diffusion model. Along one propagation direction we assume equal 
probabilities for forward and backward scattering and impose strict flux conservation in 
a stationary state situation. This leads to the following relationship between the neutrino 
density J(E) and the neutrino number flux H(E), 

H Vi { E ) 1 

^JE) -2r Vi (E) + l- ( 25) 

We can test this relationship against a complete numerical solution of the diffusion equation 
in e.g. a supernova environment where all relevant opacities are included and find agreement 
to about a factor of two. If the thermodynamical conditions and the neutrino densities 
along the propagation direction are set, relation (|A25j) defines a local neutrino number flux 
H V .[E) which in general no longer obeys flux conservation in a stationary state situation. 
Assuming that we still have a stationary state situation and that the fluxes are locally well 
represented, we can use the balance of fluxes across a infinitesimally thin layer perpendicular 
to the propagation direction to obtain an estimate of the rate R Vi of neutrinos produced in 
this layer. Denoting the propagation direction with x, we express the rate in terms of the 
prevailing neutrino density and a diffusion time scale T^ l (. with 
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R "' {E) = = KfW (A26) 

The substitution of eq. (jA25|) for H v . leads to spatial derivatives of the neutrino density 
J Ui (E) and the optical depth r Vi . As the latter is given by the negative inverse mean free 
path, -l/\ Vi (E), eq. (|MH|l can be resolved for the diffusion time scale according to 

T " (E) - — — (— &— + (^MTTyxM ) ■ (A27) 

We rewrite this estimate with a distance parameter, Ax(E), to obtain 

T%f{E) = ^S(2r„(£0 + 1), (A28) 
c 

The spatial derivative of the neutrino density in eq. (jA29|) is quite inconvenient, one would 
prefer a diffusion time scale that does not depend on neutrino densities. Moreover, the 
derivative is likely to introduce noise when evaluated in a three-dimensional numerical sim- 
ulation. Hence, we neglect this term. In physical terms this means that we assume neutrino 
sources that keep the neutrino density close to constant over a spatial interval where the 
mean free path changes significantly. This might not always be justified and is subject to 
future improvement. The expression for the distance parameter, however, greatly simplifies 
to 

Ax Ui (E) = U i (E) + l)\„ i . (A30) 



2, 

Here we recall that Ruffert et al. (1996) found the dependence 

T^f(E) = 3 AXvi{E) r Ui (E) (A31) 
c 

by calibration with a numerical neutrino transport scheme. If we go back and use r Vi ~ 1 

for the "last interaction region" to simplify eq. (jA25|) further by the approximation 

HuAE) _ 1 
cJ Vi {E) 3r Vi (EY 

we obtain eq. (jA31|) by the same analysis used to derive eq. (jA28|) . However, the distance 
parameter is then given by 

Ax Vi {E)=r Vi {E)K{E). (A32) 

In our scheme Ax defines the effective width of a layer drained by the diffusive flux, i.e. pro- 
vides the conversion between a net emitted neutrino flux (number/s/cnT2) and a production 
rate (number/s/cm"3). We choose eqs. (jA31|) and (|A32|) for our numerical simulations be- 
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cause the linear dependence in r v . allows the extraction of the energy dependence as in eq. 



(|A24jl . Approximating the neutrino distribution function in the high-density regime with a 
thermal equilibrium distribution we apply the diffusion time scale and obtain the diffusion 

rates 

^'^rwS) dE -wFkr FM <A33) 

w=rwm dE =Wk I ' FM (A34) 

with 



(Ef/) = ^K = T^M- (A35) 



(QS f ) _ T Fi{Q 

Here n Ui (E) is related to the number density by n Vi = / °° h Ui (E)dE. The statistical weights 
g Vi are 1 for v e and v e and 4 for v x . After the second equals sign in eqs. (jA33|) and (jA34|) we 
have inserted the explicit estimate (jA31|) for the diffusion time scale and (|A32|) for the dis- 
tance parameter. Note that this leakage prescription is not based on the use of mean neutrino 
energies, eq. (|A3J) is exclusively used for informative purposes. Our scheme accounts for the 
energy dependence of the neutrino opacities by integrating over the neutrino distribution. 
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Figure 1. Testing for spurious neutrino emission: Shown are the luminosities of the various neutrino species (left) and the 
corresponding mean energies. After 15 ms, corresponding to ~ 50 neutron star dynamical time scales, the total neutrino 
luminosity levels off four orders of magnitude below the emission of the full merger. 
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Figure 2. Run C (no spins, 2 x 1.4 Mq): the left panel shows the luminosities (in ergs/s) of the different neutrino flavours. 
The right panel gives the corresponding mean energies. We regard this to be the generic case. 



High Resolution Calculations of Merging Neutron Stars II: Neutrino Emission 29 




Figure 3. Run D (corotation, 2 x 1.4 Mq): the left panel shows the luminosities (in ergs/s) of the different neutrino flavours. 
The right panel gives the corresponding mean energies. 
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Figure 4. Run E (no spins, 2 x 2.0 Mq): the left panel shows the luminosities (in ergs/s) of the different neutrino flavours. 
The right panel gives the corresponding mean energies. We regard this to be an upper limit for the neutrino emission. 
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Figure 5. Shown are the SPH-particle distributions in the log(p) — T- (upper panel) and the log(p) — Y" e -plane (lower panel) 
of the generic case, run C (2 x 1.4 Mq, no initial spin; every 20th particle is displayed as a dot). For comparison with the 
supernova (SN) case we show the conditions of a collapsed 13 Mq star, 100 ms after bounce as a thick line. The particles with 
the highest luminosities are also shown: filled circles for the u e , squares for the P e and the triangles refer to the v x - 
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Figure 6. Geometry of the neutrino emission: the left column shows the total neutrino energy per time and volume in the 
orbital plane, the right column displays the vertical emission geometry. The upper two panels correspond to run C (2 x 1.4 
Mq, no spins), the middle panels to run D (2 x 1.4 Mq, corotation) and the lowest panels to run E (2 x 2.0 Mq, no spins). 
The contribution of the central object is negligible. 
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heavy mass fraction, t= 14-112 ms 
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Figure 7. Shown are azimuthally averaged values of the heavy nucleus mass fraction for matter with p > 10 10 g cm -3 (top to 
bottom run C, D and E). Although heavy nuclei are present, neutrinos from the most luminous regions (see Fig. |HJ can stream 
out vertically without encountering substantial amounts of heavy nuclei. Therefore the latter ones do not influence the total 
luminosity and the average energies. 
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Figure 8. End of run C (2 x 1.4 Mq, no spins): logarithm of the optical depths (in steps of 0.5) of the various neutrino flavours 
calculated on our grid. The first panel shows i/ e , the second P e and the third v x . The thick lines give the locus of the "neutrino 
sphere" denned as r Vi = 2/3, which is essentially the locus where the neutrinos decouple from the debris matters. For the P e 
and v x the neutrino spheres almost coincide since both are subject to scattering processes, the absorption of v e onto protons 
is only a minor correction. v e are additionally absorbed in the neutron-rich environment. Therefore the corresponding neutrino 
sphere is substantially larger. 
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Figure 10. Same as previous plot for run E. 
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Figure 11. Optical depths at a height of z= for the various neutrino flavours (top to bottom: v e , v e , v x ). 
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Figure 12. Directional dependence of the neutrino emission: the upper panel shows the neutrino luminosity per solid angle, 
(#) = A q' as a function of the angle (with the original binary rotation axis) under which the system is observed. The 
apparent luminosity is given by L^ pp = AnA Vi . The lower panel shows the rms average energies of the neutrinos as a function 
of 1?. 
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